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Abstract 

We study the generation of primordial magnetic fields, coherent over cosmo- 
logically interesting scales, by gravitational creation of charged scalar particles 
during the reheating period. We show that magnetic fields consistent with those 
detected by observation may be obtained if the particle mean life is in the 
range 10~^^ sec ^ < 10^^ sec. We apply this mechanism to minimal gauge- 
mediated supersymmetry-breaking models, in the case in which the lightest stau 
fi is the next-to-lightest supersymmetric particle. We show that, for a large 
range of phenomenologically acceptable values of the supersymmetry-breaking 
scale a/F, the generated primordial magnetic field can be strong enough to seed 
the galactic dynamo. 
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Quite homogeneous magnetic fields of intensity B ^ 3 x 10^^ Gauss are present in all 
structures of our Universe: galaxies, galaxy clusters and hydrogen clouds |l|, |], ^ One 
of the mechanisms for the generation of magnetic fields is the primordial generation of a seed 
field that is further amplified by gravitational collapse and/or dynamo [|l|. For the origin of 
the seed field, several mechanisms have been proposed recently: it has been suggested that a 
primordial field might be produced during the inflationary period if conformal invariance is 
broken Q ; in superstring- inspired models, the coupling between the electromagnetic field 
and the dilaton breaks conformal invariance and might produce the seed field ^ ^ ; gauge- 
invariant couplings between the electromagnetic field and the space-time curvature also 
break conformal invariance, but produces in general an uninterestingly small seed field 
other mechanisms are based on, for example, first order cosmological phase transitions and 



on the existence of topological defects ||T0| |TT], [12 



Recently, a new mechanism for cosmological magnetic field generation was proposed 
based on the presence during inflation of a charged, minimally coupled scalar field in 
its invariant vacuum state |T^. When the transition to radiation takes place, quantum 
creation of charged particles occurs because of the release of gravitational energy. The mean 
electric current is zero, but stochastic fluctuations around that mean give a non-vanishing 
contribution. The magnetic field induced by this stochastic current was sufficient to seed 
the galactic dynamo. However, there remained the important issue of finding a suitable 
scalar particle to generate the electric current source of the magnetic field. 

In this letter we address this problem in the context of gauge-mediated supersymmetry- 
breaking models (GMSB). In the simplest version of these models, supersymmetry-breaking 
is communicated to the visible sector through a set of massive fields, called messengers, 
which carry non-trivial quantum numbers under the gauge group |]15[ [l^ . The messengers 
$/, $7 are assumed to acquire an explicit mass Mj by the vacuum expectation value of a 
singlet field (Xf) = Mj, via a superpotential coupling 

W = (1) 

A vacuum expectation value of the auxiliary component Fj of the field Xj breaks su- 
persymmetry and induces, through the gauge interactions of the messenger fields, the 
supersymmetry-breaking masses in the observable sector. For the simplest case of sets 
of messenger fields belonging to the fundamental representation of SU(5) and a single field 
X, one gets gaugino masses 

Mi = — ^ — • (2) 

Air M ^ ^ 

where i = 1,2,3 are associated with the gauge groups U{1)y, SU{2) and SU{3)c, respec- 
tively. 

The scalar masses not affected by Yukawa couplings are given by 
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where ms are the supersymmetry-breaking masses for gauginos and scalars, respectively, 
/i is the renormahzation group scale with /i = being identified with the messenger mass 
scale, Cg is the quadratic Casimir of the scalar particle under the i-gauge group, and aj and 
bi are the corresponding gauge coupling and MSSM beta function coefficients Q. From the 
above, it is easy to see that the right-handed sleptons are the lightest scalars in the spectrum 
and, for > 1, the lighest stau can easily become lighter than the lightest neutralino. The 
lightest stau can also become lighter than the lightest neutralino due to mixing effects, for 
moderate and large values of tan/3, for any value of N. For the characteristic values of 
the supersymmetry-breaking scale F, however, the lightest supersymmetric particle is the 
gravitino. Indeed, the gravitino mass is given by 

where Mpi is the Planck scale (we are identifying F with the fundamental supersymmetry- 
breaking scale Fq). Hence, the gravitino is the lightest supersymmetric particle for any 
messenger mass M much lower than the GUT scale. 

In general, under the assumption of R-parity conservation [Q, the next-to-lightest SUSY 



particle will decay into a gravitino and a standard particle with an inverse decay rate [19 



.^iri^)V^V3xlO-GeV-, (5) 
P V mNLSP / VlOOTeV; ^ ^ 

where mNLSP is the mass of the NLSP particle and A; is a projection factor equal to the 
component in the NLSP of the superpartner of the particle the NLSP is decaying into. For 
the case of the stau decaying into a tau and a gravitino, k = 1. 

Constraints on the value of the supersymmetry-breaking scale may be obtained, for ex- 
ample, by the requirement that the gravitino density does not overdose the Universe. 
For instance, if the gravitino mass tjiq > 1 keV, the temperature at the beginning of the 
radiation-dominated epoch, called the reheat temperature T^, should be much smaller than 



the GUT scale in order to avoid overproduction of gravitinos ||2^. The exact bound on 
depends on the gravitino mass. For relatively large values of the gravitino mass, corre- 
sponding to ^/F ^ 10^ GeV, (M ~ 10^'^ GeV), an upper bound on of the order of 10^ 
GeV is obtained 0. The bound becomes even smaller for smaller values of F. On the other 
hand, for < 1 keV and for any value of the reheat temperature larger than the weak 
scale, the gravitinos will be in thermal equilibrium at early times and, for these range of 
masses, the gravitinos are sufficiently light to lead to cosmologically acceptable values of 
the relic density. 



^For more general expressions see, for instance, Ref. |Tj 

^For large values of the gravitino mass, the non-thermal production of gravitinos tends to be dominant, 
and induces a tighter bound on the reheat temperature, which may be of the order of the weak scale UM- 



3 



The reheating period can be characterized by the temperature T/ obtained by thermahza- 
tion after preheating and T^, the temperature at the beginning of the radiation dominated 



epoch [|22|. In Ref. [|T3| an inflationary model with instantaneous reheating was considered. 
In the more realistic case in which the reheating is extended in time, the number of particles 
created at the two main transitions, namely inflation-reheating and reheating-radiation, as 
well as during the reheating period itself should be calculated. 

We will work in conformal time, which is given by drj = dt/a{t). Defining r = Hrj, where H 
is the Hubble constant during inflation, and assuming that during reheating the Universe 
is matter-dominated the scale factors for the different epochs of the Universe read 

1 



inflation a/ (r) 
reheating aji (r) 



radiation (r) 
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Tj and are the temperatures of the Universe at the beginning of reheating and at the 
beginning of radiation, respectively. We have assumed that during radiation the temper- 
ature of the Universe scales with a (r) as T oc a(r)~^ while during reheating it goes as 
T oc a(r)"^ with < 6 < 1 [|4|. 

The evolution of a charged scalar field is given by the Klein-Gordon equation. If we 
expand the real and imaginary parts of the field as (27r)~^''^ / dPn 
field equation reads 



It e'- 



"+ h.c. , the 



m 



60 



a (T 



a [T 



0, 



(9) 



where k = H~^k, {k, being the comoving wavenumber) and where ^ is the coupling to 
the curvature. We will consider the mass as built up from two contributions, the zero- 
temperature mass m(0) = rrif and the thermal corrections, so that we have = m^(0) + 
gT^ (r), where g is of the order of the particle gauge coupling constants. 

For the inflationary period, we do not need the thermal corrections, as the temperature 
of that period is too low to be important. However, in supergravity theories, the possible 
presence of a non-renormalizable coupling of the inflaton field / to the scalar fields in the 
Kahler potential PT| 

Ch 



1 



3 Ml, 



(10) 



would naturally lead to a mass contribution 5m? ~ ChH^- Hence, in general, an effective 
mass of the order of the Hubble constant will be generated, although the coefficient Ch 
may be small or even zero in the case when the specific effective coupling is forbidden by 
symmetries of the theory p5 |. 
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The positive-frequency solution to Eq. (H) for inflation reads 



cPiir) = y^VT^Hi'-^ [k{l-r) 



(11) 



where H!^P are the Hankel functions, with 



16 



4 



2V 3^ 9i/2 



(12) 



where for the charecteristic values of m(0) and H during inflation, m? / H"^ ~ Ch- We will 
assume throughout this article that the scalar field couples minimally to the curvature, 
^ = 0, and that the coefficient Ch -C 1; we will briefly discuss the implications of different 
values of these quantities at the end of this article. For reheating and radiation dominance, 
we propose a WKB solution 



2\uj\ (r) 



[±i ^ u{r')dr'] 



(13) 



where 



reheating uo (r) 
radiation uj (t) 



U2 I "^(Q)^ 2( \ 



(14) 
(15) 



It is important to note that the frequency changes from imaginary to real values, at a 
certain time Tc during reheating. 

We match the solutions to the field equation in the different epochs at the transition between 
them, i.e. at the end of inflation and at the end of reheating. At both times we demand 
continuity of the corresponding modes and their first time derivatives. Care must be taken 
to match the WKB solutions through Tc, where u (r^) = 0. We obtain 



where by 



(pk (t) = ak(pt^ (r) + f3k(l)k'y (^) , 
we denote the modes during radiation and with 



(16) 
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In the above, we are ignoring the effects produced by the change of the effective Hubble con- 
stant during the inflationary period, which results in a change of the Bogoliubov coefficients 



in the far ultraviolet These effects, however, are small in the range of wavelengths 

relevant for the analysis of the generation of magnetic fields, k < ktod, where ktod is the 
comoving wave number of the relevant astrophysical scale under study (see below). 

In order to proceed with our phenomenological analysis, the values of Tj, and H in the 
previous expressions must be specified. They can be related by the age of the Universe, 
which can be well approximated by the duration of the matter-dominated epoch. This is 
given by 



ttod 



2 

3^ 



'TV 



3/2b 



T 

Tm 



AI 



3/2 



Ttod/ 



(19) 



where T^od — 10^^^ GeV is the present temperature of the Universe, Tm — 1 eV is its 
temperature at the beginning of the matter dominated epoch and for the Hubble constant 
during inflation, we shall assume that 10^^ GeV < H < 10^^ GeV. From Eq. ( [T9| ) we obtain 



Tj 



Tn(3-46)/3 

7 



3H 



rpl/2rp3/2 , 
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7 



(20) 



where the last equality stems from t 



tod 



Tic^T^ 



KPGeVV'/' 

-^7 / 



23 . For H = 10" GeV we have 



(21) 



Therefore, independently of the value of b and for the values of the cosmological parameters 
considered above, the relation T/ > is fulfilled for any value of < 10^^ GeV. 

In order to compute the magnetic field, we consider the Maxwell equation 



air 



dr 



(22) 



where a (r) is the time-dependent conductivity of the medium and j is the electric current 
generated by the charged scalar particles. Although {^j^ = 0, the two-point correlation 
function is different from zero and produces a non-vanishing magnetic field. This field 



can be expressed in terms of the two-point function of the scalar field as (see Ref. |]T3| for 
details) 



27r 



W,\„ (A) \k X kfG]^i,,^ (r„, r) G]^,,^ (r„ r') 
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cosQk (t, t') 



(23) 
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^cykm, (r) <t>t, (r') + 2a^/5,0,-^ (r) 0,^ (r') 



+ 



(25) 



with ^{t) = J'^ (It'uj {t'), {4>k-y) the positive (negative) frequency modes of the scalar 
field during radiation dominance and with (to, r) the retarded propagator for the 

electromagnetic field. W (A)^^, is the window function that filters scales smaller than A. 
Eq. (p3D then gives the magnetic energy of a field which is homogenous over volumes of 
order A^, the intensity of the field therefore being estimated as (Bj). From now on it will 
be understood A = /c^"^. where, as mentioned above, ktod is the comoving wavenumber of 
the astrophysical scale we are interested in. 

Real particle propagation can be considered as such from the moment when the frequency 
becomes real, i.e. from Tc. To evaluate Eq. (pSj) we shall proceed in the same way 
and consider only the main contribution, which originates from the last term 



13 



as m 

between brackets, which is quartic in the Bogoliubov coefficients and, within this term, 
from the non-oscillatory contributions. We perform the k integration with the same window 
function used in [^, i.e. a top-hat one. We propagate the magnetic field during reheating 
and radiation dominance until the moment of detection with the propagator given by the 
equation [d'^/dr'^ + k'^ + a (r) d/dr] Gj'fc+fc/i (tq, t) = 5 (to — r), where a (r) is the electric 
conductivity of the Universe. After the particles decay, the field propagates conformally. 
We assume that during all these periods the conductivity of the Universe is given by 



2rp 



(J r 



H 



r 



r*) 



(26) 



where a = 2b, r* = —2 and ctq = Tj/e'^H for reheating, and a = 1, = — 2 + (Tj/T^) 
and (Jo = T/^^''T^~^/^^/e^if for radiation dominance. 

For To ^ r, we have 



1/26 



Gret 
\k+k' 



(27) 



The Bogoliubov coefficients are the ones given in Eqs. (p!7| ) and (p^S]). Now we are ready 
to evaluate the time integrals in Eq. (ESI). It can be checked that the contribution from 



> 10^14 sec. We 



reheating is negligible with respect to the one from the radiation period. Also, the mass 
term dominates over the thermal correction for a particle lifetime t„ 
therefore consider the time integral 



dr' 



T^m(O) \Ti 



2H It. 



3/26' 



1/2 



(28) 



Now we are ready to evaluate the magnetic field. For this purpose it is convenient to express 
the comoving wave number in terms of the present one as ktod = i^todT-y {Tj /T^)^^^ / HTtod- 
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Replacing everything in eq. (^3]) we obtain 
Equation (EUI) gives the intensity of the field at the moment when the electric current van- 




ishes. After that, the field propagates as B^^'^if) = ^ {B'j)a {tmaxf' /O' if)'^ (i-e. magnetic 
flux conservation). Using the relation a{t) oc 1/T, the present value of the magnetic field 
is given by 



{Btn^4M)(^] (30) 

^ ^ max ' 

where T^ax is the temperature of the Universe when the particles decay, given by T^ax = 
T^a^ (r^) /a^ [jmax) = T'^^'^'^T}^~^/^^ / ^j2Htmax- Replacing everything in Eq. ( pOf ) we obtain 

lT3Phys\ _ ^ l^tod I f / \ ,3/2 /oi\ 

'- m{0)T, [tJ ^^^^ 

In the above, we have only considered the effect induced by the scalar particle and not by 
the charged particles resulting from its decay. This might seem surprising since, due to 
charge current conservation, the charged particles coming from the scalar particle decay 
might also contribute in a relevant way to the magnetic field generation. However, the 
decay of a massive scalar particle, like the stau in the case under study, will lead mostly to 
charged fermions (tau leptons, in this case) with wavelengths much shorter than the ones 
of the original scalar field. These fermions might eventually generate magnetic fields, but, 
due to the wavelengths involved, these fields will not be coherent in the scales of interest 
for our study. 

In order to apply the above formalism to the case of gauge mediated supersymmetry- 
breaking models, we should recall Eq. (^, which gives the lifetime of the NLSP, = t^ax as 
a function of the supersymmetry-breaking scale and the mass of the lightest stau. Replacing 
Eqs. (|0D and (|) into Eq. (|3l|), we obtain 

,3 7ir7/2^2 



T}^'^ [HMpif^^^ X 100 GeV 




4 

3 X lO^^GeV-^ 



3/2 

(32) 



We see that the 6-dependence has disappeared, i.e. the result does not depend on the details 
of the reheating period. Using the equivalence 1 GeV^ ~ 10^° Gauss and the numerical 
estimates H ^ 10^^ GeV, ~ 10'' GeV, Ktod — 10^'^^ GeV (for a galactic scale of the order 



of 1 Mpc), m(0) ~ 100 GeV and ^jF/k = 10^ GeV, we obtain 

{Bf y') ~ 10-1' Gauss. (33) 
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This value of the generated magnetic field is sufficient to seed the galactic dynamo, being 
also consistent with the bounds imposed by the anisotropics in the CMBR and by primordial 



nucleosynthesis 28, 29 



In the above we have given results for specific values oi y F, H and T^, for minimal coupling 
and for Ch = 0, that is for v = 3/2. It is interesting to discuss the dependence on a/F, 
Ch, as well as on departures from minimal coupling. In this case it can be checked that for 
small ktod, the Bogoliubov coefficients are given by ~ /3a: ~ 0{l)k~'^, with u given by 
Eq. (|12D. Considering a stau lifetime Tf = tmax = 10" GeV~^ and a stau mass nif ^ 100 
GeV, the value of the physical magnetic field is given by 



{Bf^') = 10(lW3)(.-3/2)-106/3+3n/2 QaUSS X 

/l(r_GeVy'~'"^/' / H y^5+8.)/12 / ^^^^ ^(5-2.) 

[ J UOii GevJ U0-38GeV. 



(34) 



Acceptable values of the magnetic field, consistent with the cosmological bounds ||2^, ^ 
are in the range 

10-^ Gauss > Bf^' > 10-^^ Gauss. (35) 

For the case u = 3/2, and for the mean values of the parameters taken above, this implies 
158/9 > n > 86/9, or, equivalently, 

10"^^ sec <Tf< 10"^ sec. (36) 

The variation of the above bounds with the cosmological parameters can be easily obtained 
from Eq. (13). 

On the other hand, for arbitrary values of u, the following bound is obtained: 

3 79 9 3 43 9 

- H n > u> - -\ n. (37) 

2 118 236 - - 2 118 236 ^ ^ 

In gauge-mediated supersymmetry-breaking models, for a stau mass rrif ~ 100 GeV and 
values of the gravitino mass rriQ ^ 1 keV, the stau lifetime, Eq. (|^), is such that 11 ^ n ^ 16, 
or equivalently 

10~^^ sec <Tf< 10"^ sec GMSB for < 1 keV. (38) 

For a mimimal stau coupling to the curvature and Ch ^ 1, that is for z/ = 3/2, Eq. (|38|) is 
in remarkable agreement with the values required to generate an acceptable magnetic field, 
Eq. di. 

The bounds on n in gauge mediated supersymmetry breaking models also imply bounds 
on u 

1.72 > z/ > 1.25. (39) 
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Comparing this expression with the value of v for minimal coupling of the scalar field, 



V ^ 3/2^1 — 4Ci^/9, we obtain that Cu < 0.68 in order to generate cosmologically relevant 
values of the magnetic field. As follows from Eq. (Q), only small modifications of the bound 
on Ch may be obtained for different values of the cosmological parameters. 

Consider now the departure from minimal coupling. Assuming that Ch ^ 1 we have 



u ~ 3/2y 1 — 16^/3. The bounds on the lifetime of the stau are satisfied for 

< ^ < 0.06; (40) 

we thus obtained for a non- negligible interval of coupling values, magnetic fields of an 
intensity sufficient for these to be cosmologically important. The upper bounds on Ch and 
^ quoted above can only be obtained for values of \/F (or equivalently n) such that the 
gravitino mass is close to 1 keV. 

The results given above were obtained for a reheat temperature ^ 10^ GeV. As we em- 
phasized above, for the range of gravitino masses we are concentrating on, the most relevant 



bound on the reheat temperature comes from Eq. (20), which assures the consistency of 
the whole approach. Larger values of the magnetic fields may be obtained by lowering the 
value of the reheat temperature. However, the final result for the magnetic field, Eq. (|3^ , 
depends very weakly on the value of the reheat temperature T^. No relevant departures 
from the obtained values would be obtained even if the reheat temperature were as low as 
~ 10^ GeV. 

In summary, we have shown that cosmologically relevant magnetic fields may be gener- 
ated by a scalar field, minimally coupled to the curvature, so far its lifetime is bounded 
by Eq. (0). The bounds on the lifetime are in excellent agreement with those obtained 
in minimal gauge mediated supersymmetry-breaking models with the lightest stau as the 
next-to-lightest supersymmetric particle, for values of the supersymmetry-breaking scale 
such that rriQ < 1 keV. This conclusion is very weakly dependent on the assumed values 
of the cosmological parameters. Moreover, contrary to many models for magnetic field 
generation proposed in the literature, the present one is related to the properties of the low 
energy effective theory and these properties can be tested in accelerator experiments in the 
near future . 
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